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ABSTRACT 

In quasi-persistent neutron star transients, long outbursts cause the neutron star crust to be heated out of 
thermal equilibrium with the rest of the star. During quiescence, the crust then cools back down. Such crustal 
cooling has been observed in two quasi-persistent sources: KS 1731-260 and MXB 1659-29. Here we present 
an additional Chandra observation of MXB 1659-29 in quiescence, which extends the baseline of monitoring 
to 6.6 yr after the end of the outburst. This new observation strongly suggests that the crust has thermally 
relaxed, with the temperature remaining consistent over 1000 days. Fitting the temperature cooling curve 
with an exponential plus constant model we determine an e-folding timescale of 465 ± 25 days, with the crust 
cooling to a constant surface temperature of kT££ = 54 ± 2 eV (assuming D = 10 kpc). From this, we infer a 
core temperature in the range (3.5-8.3) x 10 7 K (assuming D = 10 kpc), with the uncertainty due to the surface 
composition. Importantly, we tested two neutron star atmosphere models as well as a blackbody model, and 
found that the thermal relaxation time of the crust is independent of the chosen model and the assumed distance. 
Subject headings: stars: neutron — X-rays: binaries — X-rays: individual: MXB 1659-29 
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1. INTRODUCTION 

Many low-mass X-ray binaries (LMXBs) are transient, 
spending the majority of their time in a quiescent state 
with very low levels of accretion and a small fraction of 
their time in outburst where the mass accretion rate (and 
hence X-ray luminosity) increases significantly. For neutron 
stars, these repeated outbursts affect the star - the compres- 
sion of the crust due to accretion of matter induces elec- 
tron captures, neutron em issions and pycnonuclear reactions 
dHaensel & Zdunikl [1990) in the crust, which in turn heat 
the core. Over approximately 1 4 — 1 5 yr a steady-state is 
reached, in which this deep crustal heati ng during outburst is 
balanced by cooling during quiescence dBrown et al.| [T998). 

Typically outbursts in these neutron star X-ray transients 
only last weeks to months; however, in the so-called quasi- 
persistent transients, outbursts last several years. In these ob- 
jects enough heat is imparted to the neutron star crust that 
it gets heated significantly out of thermal equilibrium with 
the rest of the star, which is not the case for the major- 
ity of normal X-ray transients. Therefore, in these quasi- 
persistent sources, once the source returns to quiescence the 
crust cools down into thermal equilibri um with the cor e over 
an appreciable timescale (e.g. [Rutledge et al.l l2002bl) . So 
far, crustal cooling curves have been observ ed in two such 
sourc e s: MXB 1659-29 and KS 1731- 260 dWiinands et al.l 
1200 1L 120031 I2004t ICackett et all E>006). For both objects 
we obtained several Chandra and XMM-Newton observations 
over a pe riod of approximatel y 4 years after the end of long 
outbursts. Cackett et al. ( 2006) found that both sources cooled 
down rapidly, with the cooling well fit by an exponential de- 
cay to a constant level. This was interpreted as the neu- 
tron star crust cooling down into thermal equilibrium with 
the rest of the star. By comparison to crustal cooling models 
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of iRutledge et al.1 d2002b) the cooling curves suggest that the 
crusts both have a high thermal conductivity, and the cores 
may require enhanced levels of neutrino emission. These 
observations motivated further th eoretical study of crustal 
cooling, and models calculated by Sht ernin et al.l (120071) for 
KS 1731-260 also have a best fit with a high thermal conduc- 
tivity crust, but may not require enhanced neutrino emission 
in the core. In this paper, we present an additional Chandra 
observation of MXB 1659-29, extending the quiescent mon- 
itoring by 2.8 yr to now cover 6.6 yr after the end of the out- 
burst. In section |2] we detail the data analysis and present the 
results, while in section [3] we discuss the cooling curve and 
implications for the neutron star structure. 

2. ANALYSIS AND RESULTS 

The neutron star low-mass X-ray binary MXB 1659-29 
was in outburst for about 2.5 yr, with the source going into 
quiescence in September 2001. Since then, the source had 
bee n observed 6 times (5 with Chandra, 1 with XMM -Newton ; 
see lWijnands et al.ll2003l 2004t ICackett et al.ll2006h . On 2008 
April 27 (MJD 54583.8) a 7th quiescent observation was 
performed - we observed MXB 1659-29 for approximately 
28 ksec with the Chandra ACIS-S (ObsID 8984). The source 
was at the default aim-point on the S3 chip and was observed 
using the standard FAINT data mode. These data were re- 
duced using the latest Chandra software (CIAO version 4.0) 
and calibration databases (CALDB version 3.4.3). No back- 
ground flares were observed, so all data is used. A circular 
region centered on MXB 1659-29 with a radius of 3 arcsec 
was used to extract the source spectrum. For the background 
spectrum we used an annulus centered on the source with an 
inner radius of 7 arcsec and an outer radius of 22 arcsec. The 
source and background spectra were extracted using the CIAO 
tool psextract, with the correct rmf and arf created using 
the mkacisrmf and mkarf tools. 

MXB 1659-29 is seen to exhibit X-ray eclipses and dips, 
with an orbital period of ~7.1 hr and an eclipse length 
of 90 s dCominskv & Woodlll984l 119891: lOosterbroek et al.l 
120011) . During the eclipses we do not expect to detect 
the source and in the very first Chandra observation of 
MXB 1659-29 in quiescence, the eclipse was clearly seen in 
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th e lightcurve (|Wiinands et al.ll2003l) . Using the ephemeris 
of lOos terbroek etal.1 (l200ll) . we determined that only one 
eclipse occurred during our observation and we therefore re- 
duced the exposure time in the spectrum by 900 s to ac- 
count for the fact that we would not receive any counts during 
the eclipse. After this correction, the background-subtracted 
net count rate in the 0.5-10 keV band is (1.0 ±0.2) x 10" 3 
counts s -1 , with the corrected exposure time 26.7 ks. 

2.1. Spectral Fitting 

In this analysis we fit this latest Chandra spectrum simul- 
taneously with the previous 5 Chandra and 1 XMM-Newton 
observations of MXB 1659-29 in quiescence (we do not 
create new spectra of these previous observations here; see 
ICackett et al.ll2006j). F or spectral fitting we used xspec (ver- 
sion 1 1 . lArnaudl 1 9961) . The spectra were left unb inned due to 
the lo w number of counts and the W-statistic (Wacht er et al.1 
fl979h was used in all fits. 

Neutron star atmosphere spectra deviate from that of a sim- 
ple blackbody as they have a slightly harder tail due to the 
strong frequency depen dence of opacity for free-free absorp- 
tion dZavlin et al.ll 1 996b . This difference means that when a 
blackbody model is fitted to the X-ray spectra of quiescent 
neutron stars a significantly higher surface temperature and 
an unrealistically sm all emitting radius are determined (e.g., 
Rutledge et al. 1999). The spectra were therefore fitted with 
an absorbed neutron star atmosphere model. There are cur- 
rently a variety of neutron star atmosphere models available 
for the pure h ydrogen, low - magne tic field case that is rele- 
vant here. In [Cackett et al. (2006) we used the nsa model 
dPavlov et al.l I199U IZavlin et~ail Il996i) . which we also use 
here. In order to test the model dependence o f the cooling 
curve s we chose to also fit the nsatmos model dHeinke et alj 
2006). Finally, for completeness, we also fit a simple black- 
body model (bbodyrad). 

In the neutron star atmosphere fits, we fixed the mass and 
radius at the canonical values of 1.4 M Q and 10 km as the 
data is not of high enough quality to independently determine 
these parameters. The normalization in the nsa model is de- 
fined as 1 /D 2 , where D is the distance to the object in pc. 
In the nsatmos model there is also a distance parameter, as 
well as a separate normalization, K, which corresponds to the 
fraction of the neutron star surface tha t is emitting. We fix 
K= 1 in all our fits with nsatmos. In ICacketFe t al. (2006) 
we fixed D = 10 kpc to reduce the number of free parame- 
ters in the fit given the low number of counts in the major- 
ity of the spectra. If the distance is left as a free parameter 
the uncertainty in the parameters are dominated by the uncer- 
tainty in the distance. T he distance to MXB 1659-29 is esti- 
mated to be 10-13 kpc dMuno et al.ll200U lOosterbroek et ail 
2001) from using type-I X-ray bursts. Here we investigated 
the dependence of the cooling timescale on the assumed dis- 
tance by performing spectral fits with D = 5, 10 and 13 kpc. 
IWiinands et al.l (120041) found that when assuming 3 different 
distances to MXB 1659-29 there was little effect on the e- 
folding timescale. 

For the interstellar photoelectric absorption we used the 
phabs model, and this parameter is tied between all the ob- 
servations. Therefore in the spectral fits the free parameters 
are the absorbing column and the effective temperature. Fi- 
nally, for the one XMM-Newton observation we tie the param- 
eters between the MOS1, MOS2 and pn spectra. Results of 
the spectral fitting are given in Table Q] Note that we quote 
kT^ , the effective surface temperature for an observer at in- 



finity. The bolometric flux was estimated by extrapolating the 
unabsorbed model over the 0.01-100 keV range. All uncer- 
tainties qu£tedjindplotted in this paper are 1-a. We note 
that in ICackett et alj d2006l) the uncertainties quoted on the 
spectral fits and plotted on the figures are at the 90% level of 
confidence not l-a as is incorrectly written in the text. 

All the spectral fits show that the temperature and flux of 
this new observation has remained consistent with the last two 
observations performed approximately 1000 days previously, 
regardless of the model used, or distance assumed, in the fit- 
ting. The temperature and flux determined from the nsa fits 
with D = 10 kpc is shown in Figure Q] We found that there 
is very good agreement between the effective temperature de- 
termined by the nsa and nsatmos fits. Additionally, we 
found that the temperatures and fluxes from the neutron star 
atmosphere fits are just shifted up or down in a systematic 
way depending on the assumed distance. The temperatures 
determined from the bbodyrad fits are significantly higher 
than the neutron star atmosphere temperatures. Moreover, if 
a distance of 10 kpc is assumed, the normalization implies an 
emitting radius of ^2 km, much smaller than is realistic if the 
entire surface is emitting. Both high temperatures and small 
emitting radii are normally found w hen fitting a blackbod y 
model to quiescent neutron stars (e.g jRufledge et afl [l999). 

2.2. Cooling curves 

The results of the spectral fitting (Table Q]) show that af- 
ter the initial cooling the temperature i s consistent with be- 
ing constant over the last 1000 days. In ICackett et al.l d2006l) 
we found that the cooling curves for both MXB 1659-29 and 
KS 1731-260 could be well fit by an exponential that de- 
cays to a constant level, interpreted as the neutron star crust 
cooling back into thermal equilibrium with the core. With 
the addition of the new observation of MXB 1659-29 we fit- 
ted the cooling curve again with the same model of the form 
y(t ) = a exp[-(? - to) /b]+c. We set to to midday on the last 
day that the source was observed to be active, MJD 52159.5. 
This model was fit to the temperatures determined by the nsa 
spectral fits (for all assumed distances) as well as to the re- 
sults of the blackbody fits, allowing us to test the model de- 
pendence of the cooling curves. We do not fit curves to the 
nsatmos results as they so closely match the nsa results. 
We only fit the cooling curve to the temperatures, as from this 
fit one can directly calculate the flux via F = aT^ 4 {R a /D) 2 , 
where R a is the apparent radius for an observer at infinity. 

The results of these exponential decay plus constant fits 
are given in Table [2] and we show the temperature and flux 
cooling curves for the nsa fits (with D = 10 kpc) in Fig- 
ure Q] There is excellent agreement between the e-folding 
timescales, b, when using results from different spectral mod- 
els and assuming different distances, demonstrating that the 
cooling timescale is robust. For comparison, the cooling curve 
values for KS 173 1-260 values are a = 40 ± 3 eV, b = 305 ±47 
days, c = 70.2 ± 1.2 eV fitting the temperatures , assuming 
D = 7 kpc. These are updated from ICackett et al] d2006l) us- 
ing the correct 1-a uncertainties, and rectifying a minor error 
in the fitting code. Note, how ever, that these valu es remain 
consistent with those quoted in Cackett et al] (120061) . 

The exponential plus constant model fits the data well with 
the reduced \ 2 values all close to 1, and the % 2 probabili- 
ties, P x , all close to 0.5, as expected for a good fit. How- 
ever, the mathematical solution for the flu x from a cooling 
thin l a yer is a (possib ly broken) power-law dEichler & Chengl 
119891: [Piro et alj 120051) . not an exponential. We therefore 



Crustal cooling of MXB 1 659-29 3 



TABLE 1 
Spectral fitting parameters 



Parameter 


2688 


3794 


0153190101 


3795 


5469 


6337 


8984 




(CXO) 


(CXO) 


(XMM) 


(CXO) 


(CXO) 


(CXO) 


(CXO) 


MJD 


52197.8 


52563.2 


52712.2 


52768.9 


53560.0 


53576.7 


54583.8 


NSA,D= 10 kpc 


N H (10 21 cm" 2 ) 








2.0 ±0.2 








kT£> (eV) 


121 ±1 


85 ±1 


77 ±1 


73 ±1 


58±2 


54±3 


56±2 


fbol a (10- 14 erg cm" 2 s" 1 ) 


41 ±2 


10±1 


6.6±0.3 


5.1 ±0.3 


2.0 ±0.3 


1.5 ±0.3 


1.8±0.3 


NSA, D = 5 kpc 


Nh (10 21 cm" 2 ) 








2.9±0.1 








kT™ (eV) 


95 ±1 


69=1= 1 


63 ±1 


59 ±1 


48=1= 1 


45±2 


47±2 


F w (10- 14 erg cm- 2 s-') 


59±2 


16 ± 1 


11 =b 1 


9.0 ±0.6 


3.7 ±0.4 


2.9 ±0.5 


3.5 ±0.5 


NSA,D = 13 kpc 


Nh (10 21 cm" 2 ) 








1.6=1=0.1 








kT™ (eV) 


133 ±1 


93 ±1 


84=1=1 


79±2 


62±2 


57 ±3 


61 ±2 


F bol (10- 14 erg cm- 2 s" 1 ) 


34±1 


8.1 ±0.5 


5.4±0.2 


4.1±0.3 


1.6±0.2 


1.2±0.2 


1.5±0.2 


NSATMOS, D = 10 kpc 


N H (10 21 cm" 2 ) 








2.0 ±0.1 








*Tg (eV) 


121 ±1 


86=1= 1 


78 ±1 


73 ±1 


57±2 


54±3 


57±2 


F bol (10- 14 ergcm- 2 s-') 


39±1 


9.8 ±0.5 


6.6±0.2 


5.1 ±0.4 


2.0 ±0.3 


1.5±0.3 


1.9±0.3 


BBODYRAD 


N H (10 21 cm" 2 ) 








1.6±0.2 








Normalization b 








-i 1+0.9 
■'•-'-0.5 








kT£ (eV) 


300±11 


213±6 


193 ±5 


184 ±5 


148 ±5 


140 ±7 


147 ±5 


F bol (10- 14 erg cm- 2 s- 1 ) 


29±4 


7.2 ±0.9 


4.9 ±0.5 


4.0 ±0.5 


1.7 ±0.2 


1.3 ±0.3 


1.6±0.2 



NOTE. — All uncertainties are ler. The mass and radius are fixed in all fits to 1.4 Mq and 10 km respectively. In the nsatmos fits the normalization is fixed at 1.0. ObsIDs of each 
observation are indicated at the top, with CXO (Chandra) and XMM (XMM-Newton) denoting the observatory. 

a Bolometric flux calculated from the model over the 0.01 - 100 keV range Normalization for the blackbody model is (R/D10) 2 . R is the emitting radius (km) and D10 is (D/10 kpc) 



TABLE 2 

Cooling curve parameters from fits to the temperatures 



Parameter 


nsa 


nsa 


nsa 


bbodyrad 




D = 10 kpc 


D = 5 kpc 


D = 13 kpc 




Normalization, a (eV) 


73±2 


55 ±1 


81 ±2 


175 ±12 


e-folding time, b (days) 


465 ±25 


483 ±30 


473 ± 23 


435 ±44 


Constant, c (eV) 


54±2 


45±1 


58±2 


142 ±4 


2 

Xu 


0.8 


1.5 


1.4 


0.7 




0.52 


0.18 


0.21 


0.61 



NOTE. — The model fit to the cooling curves is of the form y(/) = aexp[-(r— lo)/b]+c, 
where fo = 52159.5. All uncertainties are lcr. 

also attempt to fit a single power-law to the data of the form 
y(t) =a(t-to) b , and find that such a model does not fit the data 
well. However, inspection of the cooling curve indicates that 
the middle section of the curve does appear to follow a power- 
law. Fitting a power-law to just these middle data points (ex- 
cluding the first and last observation) achieves a good fit. For 
fits to the temperatures from the nsa results with D = 10 kpc, 
we find the power-law index b = -0.33 ± 0.02, \\ = 0.9 and 
P x = 0.46. The effective temperature from the last observation 
is a 4.5cr deviation from an extrapolation of this best-fitting 
power-law. This again strongly indicates that the temperature 
is now remaining constant, and that the crust is thermally re- 
laxed. 

3. DISCUSSION 

We presented a new observation of the quasi-persistent neu- 
tron star X-ray transient MXB 1659-29 in quiescence, ex- 
tending the quiescent monitoring to 6.6 yr. Results from the 
first 6 observations showed that the source had cooled rapidly 
and indicated that the neutron star crust may have returned 
to thermal equilibrium with the core. This new observation 
shows that the neutron star temperature and flux remained 
consistent with the previous two Chandra observations per- 



formed approximately 1000 days before. The model depen- 
dence of the thermal relaxation timescale was investigated 
with 2 different neutron star atmosphere models as well as 
a simple blackbody model. Moreover, we assumed 3 differ- 
ent distances to MXB 1659-29. The e-folding timescales of 
the cooling curves from all the spectral fits are consistent with 
each other, demonstrating the robustness of the measurement. 
The results are cons istent with fits to the first 6 observations 
(ICackett et alJl2006h . 

With the crust thermally relaxed, we can compute the core 
temperature (here we assume D = 10 kpc). We integrate the 
thermal structure equatio n in the neutron star envelope, fol- 
lowing the calculation in iBrown et al.1 d2002l) . The inferred 
core temperature is relatively insensitive to the mass of the 
neutron star (the proper T e g increases with redshift for a fixed 
T^jj , but the surface layer becomes thinner with increasing g 
and reduces the rise in temperature in the envelope). There 
is a significant uncertainty re sulting from the de pth of the 
light element layer, however (Br own et all l2002h . We find 
the inferred core temperature to range from 3.5 x 10 7 K, for 
a deep He layer (column of 10 8 gem" 2 ) overlaying a pure 
Fe layer to 8.3 x 10 7 K for an shallow He layer (column of 
10 4 gem" 2 ) overlaying a layer of heavy rp-process ashes. We 
estimate the time-averaged mass accretion rate for the sys- 
tem to be 7 x 10" 11 M yr _1 , based on the known outburst 
behavior of 2 outbursts lasting approximately 2.5 yr, with a 
quiescent period of 21 yr. We estimated the fraction of Ed- 
dington luminosity by calculating the ratio of the average per- 
sistent oirfcuretflirxtotiie peak type-I X-ray burst flux (taken 
from lWijnands et aU 2002). From this, we estimated the time- 
averaged crust nuclear heating, assuming the heat deposited 
in the crust is 1.5 Me V per nucleon (see IBrown et al.1 fl 9981 : 
iRutledge et alJ2002ah . to be ~ 6 x 10 33 ergss" 1 , for a distance 
of 10 kpc. 

As noted previously (ICackett et al.1 120061 : iHeinke et all 
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FIG. 1. — Flux (top) and temperature (bottom) cooling curves for 
MXB 1659—29. The best-fitting model to the temperatures comprising of an 
exponential decay to a constant level is shown (solid line), where the constant 
offset is shown with a dashed line. The flux cooling curve is then calculated 
from the fit to the temperatures. The time to = 52159.5 (MJD) is the last day 
when the source was seen to be active. Error bars are lo\ Data points are 
from the nsa fits with the distance fixed at 10 kpc. 



l2007h . even for the highest core temperature compatible with 
T^, the neutrino luminos ity resulting from modifie d Urea 
cooling (for a review, see lYakovlev & P ethick 2004) would 
still be a factor of ss 30 less than the time-averaged crust nu- 
clear heating. As a check on whether there is a need for en- 
hanced cooling, we also computed the ne utrino cooling ac - 
cording to the "minimal cooling model" (Page et al. 2004), 
which includes the pair breaking and formation (PBF) neu- 
trin o emissivity, but with th e So channel suppressed follow- 
ing ISteiner & Reddyl (120081) . In this case the PBF neutrino 
luminosity from neutrons in the 3 P2 state is still sufficient to 
balance the time-averaged nuclear heating, if the core temper- 
ature is in the upper half of the range given above. Given the 



uncertainty in the depth of the light element envelope and the 
superfluid critical temperatures, we cannot exclude that the 
neutrino emission is solely due to standard cooling processes. 

The exponential plus constant cooling curve fits provide 
a measure of the thermal relaxation time of the crust. This 
relaxatio n time depends on the crust composition an d lattice 
structure dRutledge et al.ll2002bt IShternin et"ai1l2007l) . on the 
crust thickness and h ence surface gravity of the neutron star 
(Lattimer et al. 1994) , and on the distributio n of heat sources 
Shternin et al.ll2007t Horowitz et al.ll200 8b; Brown & Cum- 



ming, in preparation). IShternin et al.l (120071) showed that the 
cooling timescale in KS 1731-260 was best fit by having a 
high thermal conductivity in the crust, as if it were composed 
of a locally pure lattice. This matches molecular dynamics 
simulations (Horowitz et al. 2007, 2008a), which find that the 
dense crust plasma does indeed freeze into an ordered lat- 
tice with a high thermal conductivity. As in KS 1731-260, 
our fits to the cooling of MXB 1659-29 are again consis- 
tent w ith such an ordered, low-impurity crust. IShternin et al.1 
(2007) noted that the crust may not have completely ther- 
mally relaxed; we find a single power-law decay also fits the 
cooling curve for KS 1731-260 well, with a power-law in- 
dex = -0.12 ±0.01 and xt = 0-2 when fitting to the tempera- 
tures. Further observations are required to determine whether 
KS 1731-260 has continued to cool following a power-law 
decay or if it has rea ched a constan t kT£? indicative of a ther- 
mally relaxed crust. Sht ernin et alJ (120071) note that the cool- 
ing of KS 1731-260 can be fit without invoking enhanced 
neutrino emission in the core. If MXB 1659-29 does have 
a higher core neutrino emission than KS 1731-260 then this 
may imply that the neutron star in MXB 1659-29 is somewhat 
more massive than the one in KS 1731-260 since most equa- 
tions of state allow higher levels of neutrino emission with 
increasing central density. We note, however, that this con- 
clusion assumes that the time-averaged mass accretion rates 
in both these objects have remained in a steady state. 
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